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ABSTRACT 

Context. A new generation of powerful and efficient spectropolarimeters has recently been used to provide the first sample of magnetic 
Ap stars of accurately known ages. Modelling of these data offer the possibility of significant new insights into the physics and main 
sequence evolution of these remarkable stars. 

Aims. New spectra have been obtained with the ESPaDOnS spectropolarimeter, and are supplemented with unpolarised spectra from 
the ESO UVES, UVES-FLAMES, and HARPS spectrographs, of the very peculiar large-field magnetic Ap star HD 3 1 8 107, a member 
of the open cluster NGC 6405 and thus a star with a well-determined age. The available data provide sufficient material with which to 
re-analyse the first-order model of the magnetic field geometry and to derive chemical abundances of Si, Ti, Fe, Nd, Pr, Mg, Cr, Mn, 
O, and Ca. 

Methods. The models were obtained using ZEEMAN, a program which synthesises spectral line profiles for stars that have magnetic 
fields. The magnetic field structure was modelled with a low-order colinear multipole expansion, using coefficients derived from the 
observed variations of the field strength with rotation phase. The abundances of several elements were determined using spectral 
synthesis. After experiments with a very simple model of uniform abundance on each of three rings of equal width in co-latitude and 
symmetric about the assumed magnetic axis, we decided to model the spectra assuming uniform abundances of each element over the 
stellar surface. 

Results. The new magnetic field measurements allow us to refine the rotation period of HD 318107 to P = 9.7088 + 0.0007 days. 
Appropriate magnetic field model parameters were found that very coarsely describe the (apparently rather complex) field moment 
variations. Spectrum synthesis leads to the derivation of mean abundances for the elements Mg, Si, Ca, Ti, Cr, Fe, Nd, and Pr All 
of these elements except for Mg and Ca are strongly overabundant compared to the solar abundance ratios. There is considerable 
evidence of non-uniformity, for example in the different values of {B-) found using lines of different elements. 
Conclusions. The present data set, while limited, is nevertheless sufficient to provide a useful first-order assessment of both the 
magnetic and surface abundance properties of HD 318107, making it one of the very few magnetic Ap stars of well-known age for 
which both of these properties have been studied. 

Key words, stars magnetic fields - stars: chemically peculiar 



1. Introduction 

The existence of main sequence A- and B-type stars with strong 
magnetic fields has been known for decades. These peculiar stars 
(referred to as Ap stars) show anomalous abundances of partic- 
ular elements, for example of Cr, which may be as much as 10^ 
times (2 dex) overabundant compared to the Sun. It is also not 
uncommon to fi nd atmospheric abundances of rare earths well in 
excess of solar (Rvabc hikovai [199111 . Almost all Ap stars have 
angular momenta that are of the order of 10% or less of typi- 



* Based on observations made with the European Southern 
Observatory telescopes under ESO programmes 065.I-0644(A), 075. C- 
0234(B), 079.C-0170(A), 081.C-0034(A), and 082.C-0308(A),obtained 
from the ESO/ST-ECF Science Archive Facility. 
** Based on observations obtained at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
of Canada, the Institut National des Science de I'Univers of the Centre 
National de la Recherche Scientifique of France, and the University of 
Hawaii. 



cal values for normal stars of similar mass. The magnetic field 
strength, spectral line strengths and shapes, and brightness in 
various photometric bands all vary with the rotation period of 
the star. This variability is explained using the oblique rotator 
model: the magnetic field axis and rotation axis of the Ap star 
are not aligned with the line of sight, nor with one another, and 
several elements are distributed non-uniformly over the surface 
in a pattern that is not axisymmetric about the rotation axis. This 
means that observations through the rotation cycle lead to vary- 
ing field measurements as a result of observing the magnetic 
field from different aspects. Similarly the spectrum, and (as a 
result of line blocking and backwarming) photometric magni- 
tudes and colours, va rv as different parts of the star are observed 
( lRvabchikovalll991h . 

HD 318107 (=NGC 6405 77), a very peculiar magnetic Ap 
star, has an effective temperature T^ff = 11 800 + 500 K, lumi- 
nositv log L/Lq = 1.92 + 0.1, and mass M/Mq = 2.95 + 0.15 
(iLandstreet et al.L l2007h . It has a large global magnetic field; 
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the mean line-of-sight magnetic field component (B ) is some- 
times as large as 5 or 6 kG, and the typical size of the mean 
field modulus {B) is about 15 kG. The star is a member of the 
open cluster NGC 6405, so that it is known to have an age of 
logf - 7.80 + 0.15 (yr). This means that the st ar has completed 
about 17 + 7% of its main sequence lifetime dLandstreet et al.L 
[2007, 2008). It is still a rather young star. 

Ini t ially, a rotation period of 52.4 days was proposed by 
iNortiil (fT987h: however, more extensive photometric data ob- 
tained by Manfroid & Mathvsl (l2000l) combined with measure- 
ments of the mean surface magnetic field {B} led to the conclu- 
sion that the rotation period is f = 9.7085 ± 0.0021 days. 

The magnetic field o f this star was modelled by 
[Landstreet & Mathvs' (200(f) based on data from iMathvs et alj 
(1997), Mathys & Hubrig (1997) and furt her unpubhshed mag- 
netic field measurements by Mathys. iLandstreet & Mathvs' 
(120001) adopted a model of colinear magnetic dipole, quadrupole, 
and octopole components of strength +23700, -23600, and 
+8300 G respectively. Observations of the hemispheric ally aver- 
aged line of sight component of the magnetic field (B,) indicate 
that (B ) is always positive, which leads to the model constraint 
that the sum of the inclination of the rotation axis to the line of 
sight, /, and the angle between the magnetic field axis and the 
rotation axis, /3, must be less than about 90 degrees: i + /3 < 90°. 
This implies that the negative magnetic pole is never directly 
observed. For this model, the values for / and /? were chosen 
to be 11° and 78°, with un certainties discussed in the text of 
ILandstreet & Mathvsl (l2000l) . Note that this model does not re- 
produce the observations very well. It appears that the field struc- 
ture of HD 318107 is more complex than the axisymmetric low 
order multipole field geometry used, and is probably not even 
axisymmetric. Nevertheless, the model provides a first order ap- 
proximation to the magnetic field geometry of HD 318107. 

There are rather few hot Ap stars with strong magnetic fields 
for which detailed chemical abundance analyses have been car- 
ried out. HD 318107, with r eff = 11800 K, me asured mean 
field modulus of <B) ~ 15 kG dMathvs etal.Lll997h and v sin / « 



7 km s^ , is an excellent object for which to carry out such an 
analysis. The low v sin / of the star together with the large field 
strength make it possible to observe resolved Zeeman splitting 
in many spectral lines. New measurements of (B) and {B,} are 
reported here, with which it is possible to improve the precision 
of the rotation period and re-examine the magnetic model. 

HD 318107 is also interesting in that it has an age known 
with relatively high precision. The large sample of magnetic 
Ap stars in open clusters for wh ich magnetic fields have 
been measured and characterised by 'Bag nulo et al.l (120061) and 
[Landstreet et al. (2007, 2008) have already provided the first 
clear evidence of evolution of magnetic field strength of the fos- 
sil fields of middle main sequence stars with time (or stellar age) 
during the main sequence phase of evolution. This data set, much 
of which is based on relatively high-dispersion (polarised) spec- 
tra, also offers the possibility of studying empirically the evo- 
lution with time and stellar age of the chemical peculiarities of 
magnetic Ap stars. This is now being undertaken; HD 318107 is 
the first star from this new sample to be analysed for its chemical 
abundance and distribution properties. 

This paper will discuss our efforts to model the magnetic 
field and the chemical abundance distributions of several ele- 
ments for HD 318107. The following section discusses the spec- 
troscopic and spectropolarimetric observations used and the field 
strength values deduced from them; Sect. 3 discusses improve- 
ment of the precision of the period; Sect. 4 describes the mod- 
elling technique and the spectral line synthesis program used; 



Sect. 5 discusses the abundance models obtained; and Sect. 6 
summarizes the work presented. 

2. Observations 

2.1. New data 



Previous observations of HD 31810 7 inclu de photo- 
electric photometry repo rted by iNorthI (Il987h and by 
iManfroid & Mathvsl (l2000h and used for determination of the 
rotation period; a s eries of (unpolarised) spectra obtained by 
lMathvse tal](ll997 ') from which the mean field modulus (B) was 
measured, and used by Manfroid & Mathvs (2000) to redeter- 
mine the rotation period; and a small number of spectro polari- 
metric observations reported by iMathvs &.Hubrig'('l997) from 
which the mean longitudinal field {B^} was measured. There is 
also one new measureme nt of {B-} using ESO's FORSl, reported 
bv , Bagnulo et al.l(l2006l) . 

The spectra used in this project were obtained more re- 
cently than the observations noted above, and come from a va- 
riety of sources. The data available for spectral line synthesis 
consist of three newly acquired ESP aDOnS spe ctra and a new 
FLAMES-UVES spectrum ( Pasquini et al.L 12002 ). together with 
seven HARPS spectra and one UVES spectrum from the ESO 
archive. 

ESPaDOnS is a cross-dispersed echelle spectrolarimeter lo- 
cated at the Canada-France-Hawaii telescope, which can mea- 
sure all the Stokes parameters (/, Q, U, V). We have three flux 
(/) and circular polarisation (V) spectra which cover almost the 
entire window from 3700 A to 1.04 fim with a resolving power 
of 7? ^ 65 000. These spectra typically have signal-to-noise ra- 
tios (SNR) of the order of 200 over much of this window. Two of 
the spectra were obtained one after the other, and thus are at es- 
sentially the same rotational phase. The I spectra were used for 
measuring the mean field modulus, and for modelling the chem- 
ical element abundance distribution over the stellar surface. {B,} 
was measured from the / and V spectra and used to constrain the 
magnetic model. 

HARPS is located at the European Southern Observatory 
(ESO) La Silla 3.6 m telescope. It is a cross-dispersed echelle 
spectrograph which covers a spectral range of 3780 - 6910 A 
with R = 115 000. The available spectra have SNR between 
about 40 and 100. Since the available HARPS spectra do not 
contain any circular polarisation information, these data were 
used for measurement of (B) for improvement of the rotation 
period, and for abundance modelling of HD 318107. 

UVES is a high resolution, cross-dispersed optical spectro- 
graph located at ESO's Paranal Observatory. One spectrum was 
obtained using the FLAMES instrument, with a resolution of 
R - Al 000 in a spectral window covering the ranges 4200 to 
5160 A and 5230 to 6200 A. This spectrum was used to measure 
the mean field modulus {B). The other spectrum was obtained 
with the standard UVES red arm setting at 580 nm, using image 
slicer #3, resulting in a resolving power of 1 10000. 

The new data are summarised in Table [T] The first column 
reports the instrument used, and the second gives the Julian Date 
at the middle of the observation. The meaning of the remaining 
columns will be discussed in more detail below. 

The radial velocity of HD 318107 during each observation 
was measured either by fitting an appropriate model spectrum to 
that single observed spectrum during the process of modelling 
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Table 1. The new spectroscopic and polarimetric data used in our analysis. Recorded are the HJD and derived phases for each 
spectra. The values of {B,) for the ESPaDOnS spectra are shown followed by {B) values for all spectra measured from three spectral 
Unes: Fe ii ^6149, Nd iii /15050 and AbU5. 



Instrument 


HJD 


Phase 


SNR 






Fe II ^6149 


(B) (G) 
Nd TTT ^5050 


Nd III ^6145 

i. 1 u 111 


Window 
(A) 


RV 

(km s"') 


ESPaDOnS 


2454553.119 


0.991 


250 


4900 ± 


45 


13900 ± 380 


12600 + 340 


12400 + 280 


3690-10481 


-9 


ESPaDOnS 


2454553.150 


0.995 


250 


4856 ± 


42 








3690-10481 


-9 


ESPaDOnS 


2454652.793 


0.258 


320 


1985 ± 


31 


15400 + 380 


16300 + 340 


16800 + 280 


3690-10481 


-9 


HARPS 


2453463.898 


0.803 


80 






13600 + 420 


14900 + 410 


15100 + 420 


3783-6914 


-9 


HARPS 


2453582.645 


0.034 


90 






13600 + 500 


11900 + 680 


12000 + 670 


3783-6914 


-9 


HARPS 


2453583.693 


0.142 


40 






13900 + 420 


14100 + 550 


14200 + 570 


3783-6914 


-9 


HARPS 


2454223.642 


0.056 


70 






13200 + 310 


12300 + 510 


13000 + 410 


3783-6914 


-9 


HARPS 


2454633.613 


0.283 


55 






14700 ± 250 


16400 + 410 


16400 + 340 


3783-6914 


-8 


HARPS 


2454869.850 


0.615 


70 






13800 + 420 


15700 + 680 


15300 + 570 


3783-6914 


-8 


HARPS 


2454870.879 


0.721 


60 






13600 + 420 


14700 + 680 


15900 + 570 


3783-6914 


-8 


UVES 


2451692.894 


0.390 


220 






15900 + 210 


17000 + 280 


16400 + 380 


4795-6745 




FLAMES-UVES 


2454256.625 


0.453 


150 






16500 + 340 


16400 + 510 


16900 + 420 


4200-6200 


-8 



the chemical abundance distributions (see below), or by mea- 
suring the wavelengths of a few spectral lines of low Zeeman 
sensitivity. The measured radial velocity varies somewhat be- 
tween different methods of measuring it, as the line profiles of 
individual lines are (often strongly) distorted by Zeeman split- 
ting or by blends. However, the various methods usually differ 
by no more than about 1 km s"', so we conservatively estimate 
the uncertainty of radial velocities to be about +2 km s"'. The 
measured value for each spectrum, rounded to the nearest 1 km 
s is reported in the final column, labelled RV, of Table[T] The 
mean radial velocity of the star is -8.8 km s"', and there is 
no hint in our data, which were taken over a time span of more 
than eight years, of any variability. This result is consistent with 
the mean r adial velocity of NGC 6 405 (about -7 km s"') as 
reported bv lKharchenko et aTl ( l2005h . and supports the member- 
ship of HD 318107 in this cluster. Howeve r, our measurements 
are all quite different from that reported bv lManfroid & Mathvsl 
(l2000h who measure the radial velocity to be -1-9.4 ± 2.0 km s . 
This discrepancy may be the result of a sign error by Manfroid 
& Mathys, as our new radial velocities show no indication of 
variability. 

2.2. Longitudinal field strength measurements 

The new ESPaDOnS V and / spectra provide measurements of 
both (Bj.) and {B}. The (B^) values were obtained by using the 
first-order moment of the high signal-to-noise V line profile cre- 
ated by the technique of Least Squares Deconvolution (LSD). 
The m ethod used is discussed in detail in iLandstreet et all 
(l2008h . The line list initially used for the LSD processing was 
a generic list of about 2900 lines suitable for an Ap star of 
Tgff - 12 000 K, containing primarily lines of the iron peak ele- 
ments Ti, Cr and Fe, together with some lines of lighter elements 
such as Mg and Si. The integration limits were chosen by eye for 
each spectrum. The resulting measurements are listed in the fifth 
column of Table [T] The stated uncertainties were obtained by 
propagating errors from the underlying V and / LSD spectra. 

In addition to the three measurements repo rted here, one 
other n ew measurement of (B^) is available from lBagnulo et al] 
(I200(?). This measurement was obtained using FORSl, and sep- 
arate values are reported there for the (Z? ) value obtained us- 
ing only the Balmer lines (-1-6519 ± 55 G), and the value from 
the metallic spectrum (-1-3784 + 59 G). It is obvious that these 
two measurements agree poorly with one another. This is a phe- 



nomenon that is found in several of the measurements reported 
by Baanuloetal. (2006) (for NGC 2169-12, NGC 2244-334, 
HD 149277, CD -48 11051, and HD 318107). In each case the 
star has a (B,) field of several kG, and all have T^ff values above 
about 1 1 500 K. Because the field values are very large, there 
is no question about the validity of the detections, but the {B^} 
values measured in different ways can differ by as much as a 
factor of two. Why this kind of discrepancy occurs for some (but 
not all) large-field stars is not yet understood. It might be con- 
nected with the failure of the weak-field approximation which is 
used to deduce (B,) values from FORS spectra. However, it may 
simply be the consequence of magnetic field vectors which vary 
greatly over the visible hemisphere, possibly even on a fairly 
small scale, and which are sampled differently when averaged 
using different chemical elements, which in general are not dis- 
tributed uniformly over the stellar surface. 

To explore the possibility that inconsistent field measure- 
ments are due mainly to the inhomogeneous structure of the 
magnetic field and the surface distribution of various elements, 
we have remeasured the longitudinal field (B,) in our three 
ESPaDOnS spectra using spectral lines of only a single element 
at a time. There are enough lines of Si, Ti, Cr and Fe to provide 
LSD field measurements of high enough accuracy to distinguish 
differences at the few hundred G level. The results of this exper- 
iment are fisted in Table|2] 

Table 2. Measurements of {B^} using LSD masks composed 
of lines of single elements (Si, Ti, Cr, and Fe) from the three 
ESPaDOnS spectra (cf Tableffl). 



JD (S.XSi) <S,>(Ti) <S,>(Cr) <S,>(Fe) 

-2450000 (G) (G) (G) (G) 

4553.119 7031 +379 4193 + 231 4215 + 139 5430 + 50 

4553.150 7107 + 370 4222 + 221 4124+ 132 5237 + 47 

4652.793 2007 + 269 708 + 148 878 + 92 2417 ± 32 



From Table |2] a number of conclusions may be drawn. First, 
the very satisfactory agreement of the individual entries in the 
top two lines, derived from two independent spectra obtained 
within an hour of one another indicates that the claimed pre- 
cision (measured by repeatability) of the field measurements is 
believable. On the other hand, in each of the ESPaDOnS spec- 
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tra, the values of {B,} measured using different elements reveal 
remarkable differences. Cr and Ti seem to sample the field in 
similar ways; both near phase and phase 0.25 the measured 
(Bj.) values are not significantly different. Fe, however, seems to 
sample more strongly regions of field lines parallel to the line of 
sight than the other two iron peak elements, and this tendency 
is still stronger in the sampling of the field associated with Si, 
which, however, does not sample the field in the same way as 
Fe. 

The measurements of the field using different elements thus 
reveal exactly the same general phenomenon already found in 
the comparison of the {BJ value obtained using lines of H to 
the values obtained using metal (primarily Fe) lines in the sin- 
gle field measurement with FORS 1 . These results strongly sug- 
gest inhomogeneous distributions of various elements, perhaps 
combined with significant field structure variation on a small or 
medium scale, have large enough amplitude to greatly affect field 
measurements. 

2.3. Field modulus measurements 

The ESPaDOnS, HARPS and FLAMES-UVES spectra all show 
spectral lines with resolved Zeeman splitting. We have made 
measurements of the mean field modulus (B), which we have 
compared to those of Manfroid & Mathys (2000) in order to im- 
prove the precision of the rotation period of HD 318107. The 
measurements were made using the spUtting observed in the Fe ii 
.16149 line, and in Nd iii ^5050 a nd /16145. The Fe ii . 1 6149 a nd 
Nd III /16145 lines were used by Manfroid & Mathvsl (l2000t) to 
measure {B} and so to permit direct comparison are used in this 
study as well. In addition to these two lines, we also measured 
the Nd III /15050 line because its Zeeman pattern is very simi- 
lar that of Nd m /16145. In the stellar spectra this line has three 
clearly distinct groups of Zeeman components, and because the 
Lande splitting factors for the lower and upper levels are almost 
equal, the actual detailed splitting pattern does not deviate much 
from a simple triplet. Thus Nd iii /15050 is an excellent line for 
measuring (B) in order to have a check on our measurements 
of the /16145 line. We were able to take advantage of this line 
because of the the large spectral windows covered by our new 
spectra. 

The atomic data required to convert splitting measure- 
ments into (B) values were taken from t he Vienna Atomic 
Line Database (VALD^ dKupka et al.L Il999t iRvabchikova et al.L 
[19971 iPiskunov et all 1 1 995h . except for some lines of Si, which 
were given Lande factors assuming L-S coupling. To obtain {B} 
values, the equation 

AA 

^^^%.67xl0-i3.l2z' 

was used, where AA is the shift of the cr component from the zero 
field wavelength, Aq is the location of the spectral line with zero 
field, and z is the effective Lande factor of the line (respectively 
1.17, 1.00, and 1.35 for .15 050, 6145, an d 6149 ). For a more 
comp lete discussion, refer to lMathvs et alj (Il997h or lLandstreetl 
(1200 9). The goal was to reproduce the methods used by Mathys 
and collaborators as closely as possible so as to have (B) values 
that could be directly compared to theirs, in order to use the new 
measurements to improve the precision of the rotation period. 
For the Fe ii line, the two sigma components of the Zeeman pat- 
tern are separated by as much as 1 A and are easily measurable. 
The Nd iii lines have three components and the shifts between the 
central and blue-shifted components were measured to derive the 



mean field modulus. The sigma components of both Nd iii lines 
are nearly symmetric, which allows for the measurement of (B) 
to be nearly approxima tion free. This is c onsistent with the pro- 
cedures outlined bv Manfroid & Mathvsl (j2000). The results of 
these measurements are summarized in Table [T] in columns 6 
through 8. 

To estimate uncertainties, multiple measurements for each 
spectrum were made by repeatedly fitting Gaussians to indi- 
vidual Zeeman components of measured lines, using the fit- 
ting capability of the IRAF splot function. The location of the 
Zeeman components that determine AA can be determined to 
within +0.01 A or less. The corresponding AA was then cal- 
culated for each measured line. The measurements made from 
Equation ([TJ are extremely sensitive, with deviations in AA of the 
order of 0.01 A changing the measured (B) values by as much 
as 800 G (depending on the spectral line being measured). 

It is reassuring that the measurements of {B} for both Nd iii 
lines agree within the estimated uncertainties; again, this sug- 
gests that the uncertainties are realistic. We notice in the (B) 
measurements the same general phenomenon observed in the 
{B,} measurements, namely that the field strengths as measured 
from a single spectrum using lines of different elements are not 
the same. The difference is not as striking for (B) as for {B.}, but 
presumably points to the same phenomenon: the abundance dis- 
tributions of the elements used, and perhaps the field structure, 
are sufficiently inhomogeneous on a scale smaller than a stel- 
lar radius that different elements sample the field geometry in 
significantly different ways. However, it is also remarkable that 
the value of (B) varies through the rotation cycle by only about 
+ 10 - 15% from the median value observed. 



3. Refinement of the rotation period 

iManfroid & Mathvsl (l2000h used photometric data, and (B) mea- 
surements taken between August 1992 and September 1998, 
to test the rotation period of HD 318107 proposed by iNorthI 
(Il987h . They found that the the 52.4 d period did not fit all 
their data. They showed that instead a period of P = 9.7085 + 
0.0021 d fit both the photometry and (B) measurements. When 
we use this period, we choose the zero point for phase to be 
JD 2449397.828 wh i ch is equivalent to the zero point used by 
Ma nfroid & Mathvsl (l2000l) . but close to the mid-point of their 
{B} dataset. We find that the phases of our recent data have a 
phase which is uncertain by approximately 0. 1 cycle relative to 
the older measurements (this is the phase uncertainty computed 
from the one-sigma period uncertainty). Our new {B~} and (B) 
data, and the / spectra we wish to model for abundances, can- 
not be accurately phased with the older measurements without 
improving the precision of the period. 

In principle, the period can be improved both by fitting 
our new (B ) measurements to the five measurements used by 
Landstreet & Mathvs (2000), and by optimising the fit of our 
new (B) data to the well-determined variations of this field mo- 
ment as measured by both the Fe n and the Nd iii line. In practice, 
however, the variation of {B^} with phase is too sparsely sampled 
by both the older field measurements and our own for compari- 
son to allow us to refine the period. This problem is made more 
difficult for HD 318107 by the very different values of that 
are found at any particular phase by only moderately different 
measurement methods. 

Hence we turn to the new (B) data. Fortunately the numerous 
(B) values as measured by Fe ii 6149 A and by Nd m 6 145 A vary 
in significantly different ways, and from IManfroid & Mathvi 
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(I2OOOI) we have convenient polynomials describing the varia- 
tions of (B ) with phase for these two elements. Note that in 
Table 2 of Manfroid & Mathvsl (l2000l) the coefficients A of 
the polynomials for the two spectral lines have been inadver- 
tently exchanged (private communication with Mathys). Varying 
the period through the range allowed by the uncertainty of 
+0.0021 d, we find that our new (B) measurements fall very 
nicely on the appropriate polynomials for a limited range of 
periods. The best fit as judged by the of the fit for Fe 11 
is P = 9.7087 + 0.0012 d, while both Nd m lines agree on 
P - 9.7089 + 0.0007 d, for a conservative global average of 
P = 9.7088 + 0.0007 d, and an adopted zero point of JD 
2449397.828. This reduces the phase uncertainty of our new data 
with respect to the data from the 1990s to about +0.03 cycle. All 
phases reported here are on this system. 



A comparison of our n ew (B) data to the polynomials of 
iManfroid & Mathvsl (l2000h is shown in the upper two panels of 
FigurelU It is clear from this figure that the agreement of our data 
with the previous {B) measurements is very satisfactory, and it is 
also obvious why significant phase shifts of especially the Nd iii 
data rapidly lead to decreased quality of the fit between the two 
data sets, so that the new period is accurately determined. 

The improved period is then used to phase the older (B,) data 
shown by Landstreet & Mathvs (2000) with the new measure- 
ments. The results are shown in the bottom panel of Figure [1] 
Several points emerge from this figure. First, the new data are 
generally consistent with the older measurements provided that 
we use the metal line (B,) value from FORSl at phase 0.91. 
At this phase the Balmer line field strength is about 2500 G 
higher than the older measurement at almost the same phase. 
Secondly, the variations of both (B) and (B.) (from <p - 0.5 
to (j) = 1.0) appear to be rather non-sinusoidal. These facts, as 
discussed above, very strongly suggest that either the detailed 
metal abundance distributions or the field geometry, or probably 
both, are complex. This star appe ars to belong to the small sam - 
ple of stars, including HD 32633 (Gl agolevskii & Ge rth', '2008[). 
HD 37776 (Kochukhov et al., 2010), HD 133880 (Landstreet, 
1990), HD 137509 (lKochukhovt ,2006). and r Sco (Don ati etal . 
2006) in which the magnetic field may depart in an important 
way from the generally dipolar topology of the magnetic fields 
of most Ap stars. (The field and surface He distribution of the 
outstandingly peculiar magn etic He- strong star HD 37776 have 
recently been mapped by Kochukhov et alj (|2010|) . and the re- 
sulting map reveals a quite astonishing degree of complexity.) 
Thirdly, the fact that the FORSl metal line {B,} value agrees 
with the older metal line measurement, while the Balmer line 
(B,) measurement does not, is certainly consistent with the idea 
that the difference between the two FORS 1 (B ) values is due to 
the combination of a rather complex field distribution and sub- 
stantial differences in the way that the field is sampled by the pre- 
sumably uniformly distributed H and the strongly patchy metal 
distributions. It is clear that it would be extremely interesting to 
obtain V (and probably Q and U spectra) of this star with rather 
fine phase resolution, at least of order 0.05 cycle, to make it pos- 
sible to explore more fully the field and abundance distribution 
structure. 

4. Modelling technique 

To model the spectrum of HD 318107 we use the rn a gnetic 
spectrum synthesis programme ZEEMAN dLandstreeti Il988t 



iLandstreet et"an. Il989h . This programme is designed to com- 
pute the emergent spectrum of a star of specified T^ff and log g, 
with a specified magnetic field strength and geometry (currently 
characterized as a sum of colinear dipole, quadrupole, and oc- 
tupole, at a specified angle /3 to the rotation axis). Either a uni- 
form surface abundance distribution or a distribution which is a 
simple function only of latitude in the frame of reference of the 
magnetic axis can be specified. The parameters of the magnetic 
field model are chosen to match the observed field strength mea- 
surements o f (B^) and (B) as a function of phase, generally as 
described bv ILandstreet & Mathvsl (l2000l) . The magnetic model 
can be further tested by comparison of computed line profiles 
to observed ones in cases such as HD 318107 where line split- 
ting is an important line broadening mechanism. In principle, 
since ZEEMAN computes all four Stokes parameters, compari- 
son could be made with observed polarised spectra, but this pos- 
sibility is not implemented at present. 

The abundance variation with magnetic latitude is currently 
specified in the form of uniform abundance values on each of 
one to six rings of equal extent in latitude. The abundance of one 
element at a time, and its variation with latitude, can be varied by 
the programme to optimise the fit to a set of spectral line profiles 
in one or several spectra taken at various rotational phases. The 
output is then a coarse model of the abundance distribution of the 
element in magnetic latitude, effectively a very low resolution 
ID map. This modelling is repeated for various elements. 

It is important to understand that the abundance model used 
by ZEEMAN was adopted to describe the kind of abundance 
variations observed in 53 Cam (HD 65339), a magnetic Ap star 
in which there are huge differences in large-scale mean abun- 
dances bet ween the two hem ispheres surrounding the two mag- 
netic poles (iLandstreetL 11988 ). In this case, a first order descrip- 
tion that accounts only for abundance variation with magnetic 
co-latitude is an appropriate approximation. As will be discussed 
below, such large-scale variations does not seem to be a promi- 
nent feature of HD 318107, so this particular model is not a 
good approximation to the true abundance map. However, us- 
ing it does appear to provide some useful information about the 
nature of the abundance distributions. 

ZEEMAN includes a reasonably complete model of LTE line 
formation and radiative transfer in a magnetic field. The stellar 
atmosphere model is interpolated from a grid of precomputed 
solar abundance ATLAS 9 models. The four equations of radia- 
tive transfer for polarised light are solved on a grid with 0.01 A 
spacing. Local line profiles are computed as Voigt profiles based 
on atomic parameters taken from the VALD database or rea- 
sonable approximations. Line blending is correctly computed, 
by adding (polarised) line opacities due to various lines before 
solving the equations of transfer The correctness of many of the 
oscillator strengths used has been confirmed by synthesis of a 
large number of non-magnetic stars (for which the abundance 
of each element is essentially constant over the surface). Note 
that the modelling uses atmosphere models computed for solar 
abundance throughout; the atmosphere model is not updated as 
the abundance analysis proceeds. With rather large overabun- 
dances of some Fe peak elements found in this analysis, this fact 
certainly compromises the accuracy of the resulting abundance 
values (Khan & Shulyak, 2007). 

The models of the field geometry (a low order axisymmet- 
ric multipole expansion) and of the distribution of elements over 
the surface (a few simple discrete rings axisymmetric about the 
magnetic axis ) are quite schematic. Co mpared to the sort of map- 
ping done bv , Kochukhov et al.1 (I2004r2()10 ). who map field vec- 
tor and element distributions over the surface with ~ 20° spatial 



6 



J. D. Bailey et al.: Magnetic Field and Abundance of HD 318107 



03 



17000 
16000 
15000 
14000 
13000 



V 



FeII-6149 
MM 2000 



03 



18000 
17000 
16000 
15000 
14000 
13000 
12000 
11000 



Til 



ri- 



. T 5 / 

■t-r 



NdIII-6145 ^— H- 
NdIII-5050 I — 



6000 



5000 



4000 



o 

m 3000 



2000 



1000 



i \ 



V-' 



^1 



-0.2 0.2 0.4 0.6 0.1 

Phase 



1 1.2 



Fig. 1. The t op two plots show our new {B) data phased with our best fit period, superimposed on the polynomial fits of 
iManfroid & M athvs (2000) to their larger {B) dataset (green dashed curve). In the top panel the second, slightly different, {B) curve 
(red dash-dotted line) is the variation of {B) predicted by our adopted field ge ometry. The t op panel data shows our individual mea- 
surements of {B) using Fe ii A6\A9 (data with error bars are ours, dots are from lManfroid & Mathvs (2000)); the second panel shows 
our measurements of {B) using Nd iii /15050 (filled symbols) and /16145 (open symbols). The bottom figure depicts the {B ~) field 
variati ons observed in HD 318107. The data points in green (open squares) are the measurements shown by L andstreet & Mathvsl 
(l2000h . shifted in phase to the pr esent system; the red p oints (filled squares) are the new data from Table [T]plus (at phase 0.91) the 
metal line (B.) value reported bv lBagnulo et al.l ( l2006h . The dash-dotted line denotes the (B.) rotational variation predicted by our 
adopted field geometry and phase system. 



resolution in both surface coordinates, our models, in general, 
represent only a very rough approximation. However, the data 
set required for detailed mapping is far more extensive than the 
few, irregularly distributed, mostly / spectra, obtained with a va- 
riety of resolving powers, and covering a large range in SNR, 
that are available at present for HD 318107. The type of mod- 



elling we carry out here is quite appropriate for the limited data 
available, and provides a first exploration of physical conditions 
and chemical abundances in the atmosphere of this star, which 
can be used for statistical studies in which each star is charac- 
terised by only a few numbers, and for identifying stars of such 
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great interest that they warrant the extensive investment in ob- 
servations needed for detailed mapping. 

4.1. Magnetic field model 

Before starting abundance analysis, we need to establish an ap- 
proximate magnetic field model to use. An important quality of 
this model should be to reproduce reasonably well the phase 
variations of Zeeman splitting, so that the computed spectral line 
profiles are similar in form to the observed lines, and model fits 
are primarily sensitive to chemical abundance. 

To estimate the inclination / of the rotation axis and the pa- 
rameters ofa suitable mode l magnetic field structure, we follow 
lPrestonl(ll967l[T9"7Q) and .Landstreet & Mathvsl (l2000t) . 

An estimate of the radius of the cluster member HD 318107 
can be made from the valu es of logCTpff) = 4.072 + 0.02 and 
log(L/Lo) = 1 .92 + 0. 1 from lLandstreet etal.1 dlOOT). The stellar 
radius is then found to be /J - (2.22 ± 0.34)7Jo. Our adopted 
rotation period of P = 9.7088+0.0007 days (Sect. 4) and v sin / = 
7 + 2 km s"' (see below) may then be used in the equation 



where vsin/ is in km s and R is in solar units, to obtain / - 
37° + 15°. This is not a very strong constraint, but provides a 
useful check to the value of / = 22° derived below. 

From the form of the (Br) and (B) variations with rotational 
phase, we can get a rough idea of the parameters the multipole 
model of the field will require. The fact that both magnetic mo- 
ments (especially {B.}) vary significantly with phase indicates 
that i must be substantially different from zero, as we found 
above. The fact that all the {B,} data have the same sign indi- 
cates that, as the line of sight executes a cone around the rotation 
axis on the star, it always remains within one magnetic hemi- 
sphere, so / +yS < 90°. We expect that the sum of the polar fields 
will be of order 15 kG. From the phase behaviour of {B^}, it ap- 
pears that the closest approach of the line of sight to the magnetic 
axis occurs at about phase 0.9, even though the value of (B) is 
near minimum there; apparently the local field strength increases 
fairly rapidly into the magnetic hemisphere that is observed only 
briefly around phase 0.4. 

To obtai n a d etailed model the FORTRAN programme 
FLDSRCH (iLand street & Mathys, 2000) was used. The pro- 
gramme takes as input values of (B;) and (B) at four equidistant 
phases. The relative importance of these two kinds of data can 
be varied by adjustable weights. Since we are particularly con- 
cerned with having a field model that predicts the correct amount 
of Zeeman splitting in lines at various phases, we have given a 
relatively large weight to the (B) data. FLDSRCH calculates (Z? ) 
and (B) as functions of phase for a field composed of colinear 
dipole (Bd), quadrupole (Bq), and octopole (Boct) components. 
The programme iteratively searches for the values of /, /3, B^, 
Bq, and Boct that best fit the values (B-.) and (B) provided as in- 
put. Along with the best fit parameters, the programme provides 
calculated values of (Bj> and (B) as functions of phase. 

The best fit model found in this way (with of course some 
substantial arbitrariness because the data are clearly more com- 
plex than can be reproduced by the simple magnetic field model 
adopted) are / - 22°, /3 - 65°, and polar field strengths of 
Bd = +25600 G, Bq = -12800 G, and Boet = +900 G. The 
{B-){(p) and {B){(p) predicted by this model are plotted in Figure 
[1] and are in reasonably good agreement with observations of 
(B> but reproduce only the general amplitude of the (Bj) vari- 
ations. This model is a poor representation to the actual field 



structure, thus specifying uncertainties would give an unrealistic 
impression of how accurate this model is. However, experimen- 
tation with fits to the variations of {B-) and (B) suggest that we 
can best reproduce the observed data when the model parameters 
are in this region of parameter space. 

The resulting geometry is shown in Figure |2] The figure is 
a graphical representation of the oblique rotator model that dis- 
plays the relationship of the line of sight to the two angles / and 
P and to the field axis for the phases 0.4 (left figure) and 0.9 
(right), when the line of sight is respectively farthest from and 
nearest to the visible magnetic pole. 



4.2. Chemical abundance model 

With the stellar parameters and a model of the magnetic 
field established, the elemental abundance distributions are ob- 
tained using the synthesis program ZEEMAN described above. 
ZEEMAN compares the calculated line profiles to the observed 
profiles and varies the abundance to minimize the reduced 
of the difference between the model and the observed / spectral 
lines. The programme automatically adjusts v sin / and the radial 
velocity vr to optimal values at each iteration. For the analy- 
sis here, T^s was set to 1 1 800 K, and log g to 4.2, consistent 
with the value 4.22 + 0. 13 determined from M - 2.95 ± O.15M0 
dLandstreet et all l2007h and the radius found above. Note that 
because only one ionisation stage is detected for almost all ele- 
ments, and because the shapes of spectral lines are fit only ap- 
proximately at best, we do not obtain very strong constraints 
on the value of T^s from modelling the observed spectrum. 
However, w e note that our adopted Teff is consistent with the 
model fits o f lKuruczl(ll979h when comparing the computed and 
observed H/? profiles. 

Atomic data were taken from the VALD. The best-fitting 
value of vsin/ was found to be 7 + 2 km s"'. The best fit re- 
quired some experimentation and visual evaluation, as the v sin / 
broadening is smaller than the Zeeman splitting, so the optimal 
value is best found by looking at how various v sin / values blend 
the Zeeman components together. The microturbulence parame- 
ter was set to km s"', as the strong magnetic field is expected 
to completely suppress the convective motions underlying this 
line broadening agent. 

With the twelve spectra listed in Table[T] we have fairly good 
phase coverage of the rotation cycle of HD 318107. The largest 
gap is between phases = 0.80 and 0.99. However, some of the 
HARPS spectra have quite low SNR. Considering that none of 
the spectral lines seem to vary very much in equivalent width, 
it appears that an adequate sampling of the rotation cycle for 
this first exploratory modelling would be to start with a consid- 
erably smaller number of spectra, focussing on the best avail- 
able ones. We have chosen to model only five spectra: the two 
ESPaDOnS spectra at = 0.995 and 0.258, the UVES spectrum 
at (/I = 0.390, and two HARPS spectra at <^ = 0.615 and 0.803. 
For lines below 4795 A, we replace the UVES spectrum with the 
lower resolution FLAMES-UVES specti'um (cj) = 0.453). With 
this choice of spectra to model, we are using the spectra with the 
best resolving power and SNR, we have roughly uniform spac- 
ing of about A(p ~ 0.2 in phase (and therefore sample coarsely 
the full range of available directions of the line of sight relative 
to the magnetic axis), and we reduce the computation time con- 
siderably. 

Spectral lines respond non-linearly to changes in abun- 
dance. Thus fitting several lines of different strengths provides 
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Fig. 2. The geometry adopted for HD 3 1 8 1 07 . The vertical axis is the rotation axis of the star. The angle between the line of sight and 
the rotation axis is / = 22°. The angle between the rotation axis and the magnetic field axis is /? = 65°. For HD 318107, i + p < 90° 
and therefore, as the star rotates, the negative magnetic pole is never directly observed. The left panel depicts an orientation where 
the magnetic axis is farthest from the Une of sight {<f> = 0.4) and the right panel is when the magnetic axis is nearer aligrmient with 
the hne of sight ((f> = 0.9). The bands used for dividing the abundance distribution are shown as solid Unes and are normal to the 
magnetic field axis. 



a stronger test of the correctness of the model than fitting one 
single Une. Where possible, the best-fit abundance distributions 
were found by simultaneously fitting multiple unblended lines 
of a given element within a window of a few 10s of Angstroms. 
The spectral lines modelled are discussed in the sections on in- 
dividual elements. 

HD 318107 does not fulfill the conditions for which Zeeman 
was designed (see Sect. 4). The line of sight remains firmly in a 
single hemisphere, and in fact, because the / value of our model 
is small, the Une of sight only varies in magnetic co-latitude over 
a total range of about 45°. Furthermore, the strength of most 
spectral lines vary with rotational phase only rather modestly, 
suggesting that the pole-to-pole variations of abundance may not 
be very large. 

Nevertheless, we experimented extensively with multi-ring 
models (i.e. we allowed the abundance to vary with magnetic 
latitude when we searched for the best fit abundance model for 
each element). We found that a maximum of three rings should 
be used (abundances constant on rings 60° wide measured from 
the visible magnetic pole) in order to achieve convergence to 
the same model regardless of the initial abundances assumed. 
When we ran 3-ring models using different spectral regions (and 
thus diflferent sets of spectral lines), we found for a few elements 
(particularly Si and Ti) that the best fit was for a model with 
substantial differences (~ 1 dex) between one ring and another, 



but no consistent pattern emerged from modelUng diflferent spec- 
tral regions. We eventually concluded that the best fit models 
with abundances which seem to vary strongly with magnetic co- 
latitude were probably actually fitting some other kind of abun- 
dance inhomogeneity, probably on a smaller scale, as already 
strongly suggested by the very different values of {B^} found at 
a given phase using Unes of difl'erent elements (cf. Table 2). 

Our conclusion was thus that the models found by ZEEMAN 
with strong variation in abundance with magnetic co-latitude 
are unphysical. In fact there is no convincing evidence for 
strong abundance variations on a global scale. We have there- 
fore adopted a much simpler approach to abundance modeUing. 
The abundance of each element studied is assumed to be uni- 
form over the stellar surface. The best fits we obtain to groups 
of Unes in various spectral windows of our chosen five spectra, 
sampUng the rotation cycle of the star, of course do not yield 
the same mean abundance in aU windows. We use the standard 
deviation of the abundances found in different windows as a sim- 
ple measure of the uncertainty of the value of the best fit mean 
abundance. It will be seen that the uniform abundance values de- 
termined for different spectral windows are roughly concordant, 
and that with this simple assumption we can fit all spectral lines 
modelled fairly well at all phases, giving us some confidence that 
the derived average abundances are meaningful. 
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Fig. 3. Spectrum synthesis of the region 4795 - 4815 Ausing mean abundances. Black curves are observed spectra, red curves are 
models. From top to bottom phases 0.991, 0.258 (both ESPaDOnS), 0.390 (UVES), and 0.615 and 0.803 (both HARPS). The line 
of sight is closest to the magnetic pole at phase 0.9, and to the magnetic equator at phase 0.4. Note that some strong Unes are still 
unidentified. 



5. Chemical abundances of individual elements 



Table 3. Mean abundances of elements studied 



Because of the large wavelength coverage of the available 
spectra of HD 318107, our data are sufficient to obtain first ap- 
proximations to the chemical abundance distribution for O, Mg, 
Si, Ca, Ti, Cr, Fe, Nd, and Pr, and useful upper limits for the 
abundance of He and Mn. The large wavelength coverage also 
means that multiple lines of many elements can be found that 
have a variety of strengths and splitting patterns. 

The mean abundance values found for the elements stud- 
ied are listed in Table |3] together with the uncertainties esti- 
mated from the scatter of fits to (usually two) different spectral 
windows. For elements for which only a single line or region 
was available, we have estimated the uncertaint y as ±0.2 dex. 
For comparison, the solar abundances reported bv lAsplund et al.l 
(1200 9) are tabulated as well. The quality of the fits is illustrated 
with short spectral windows, showing the comparison of uniform 
abundance models with the five observed / spectra used for the 
fitting, in Figures[3]and|4] In each figure the model spectra have 
been computed with the global mean from Table |3] rather than 
the local best fits. The lines used and the consistency of the fits 
are discussed element by element below. 

5.1. Helium 

The strongest lines of He i in the available spectra should be 
the triplet lines at 4471 and 5876 A. Neither of these lines is 
unambiguously detected, and fitting these regions provides an 



Element 


log(Mx/WH) ± cr 


\og(nx/niih 


He 


< -2.50 


-1.07 


O 


-3.65 + 0.2 


-3.31 


Mg 


-5.70 ±0.2 


-4.40 


Si 


-3.65 ±0.15 


-4.49 


Ca 


-5.55 ± 0.2 


-5.66 


Ti 


-5.05 ±0.1 


-7.05 


Cr 


-4.50 ±0.15 


-6.36 


Mn 


< -6.00 


-6.57 


Fe 


-3.00 ± 0.2 


-4.50 


Pr 


-6.35 ± 0.2 


-11.28 


Nd 


-6.20 ±0.15 


-10.58 



upper limit to the abundance of helium of at least a factor of 20 
below the solar abundance. HD 318107 is definitely a He-weak 
star 



5.2. Oxygen 

The nine weak O i lines at 6155-56-58 A are blended with lines 
of Pr, and are too weak to model usefully. In contrast, the intrin- 
sically stronger lines at 7771-74-75 A are clearly visible and ap- 
pear to be unblended. However, this spectral region is only avail- 
able in the ESPaDOnS spectra, and so we have determined the 
mean abundance using only two spectra. The blend of Zeeman- 
split features is rather well fit in the two spectra about 0.25 cycle 
apart, and apart from possible non-LTE effects in these lines, the 
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Fig. 4. Spectrum synthesis of the region 5038 - 5058 Ausing mean abundances. Black curves are observed spectra, red curves are 
models. From top to bottom phases 0.991, 0.258 (both ESPaDOnS), 0.390 (UVES), and 0.615 and 0.803 (both HARPS). The line 
of sight is closest to the magnetic pole at phase 0.9, and to the magnetic equator at phase 0.4. 



O abundance is unambiguously determined, and appears to be 
about 0.3 dex lower than the solar value. 



5.3. Magnesium 

The only clean line found of sufficient strength to model is Mg n 
/1448 1 . This line is reasonably well fit at all phases with a uni- 
form abundance, and indicates that Mg is about a factor of 20 
lower in abundance than in the Sun. 



5.4. Silicon 

In the available spectra, the cleanest reasonably strong lines are 
the single Hne at 5041 A, the doublet at 5055-56 A, and the sin- 
gle lines at 5955 and 5978 A. The fit to the two longer wave- 
length lines is quite good at all phases with a single mean abun- 
dance. The fit to the shorter wavelength lines is a little less sat- 
isfactory, as will be seen in Figure |4] The fit to these lines is 
roughly equally good at all phases, suggesting that no large-scale 
abundance variation from hemisphere to hemisphere is detected. 

5.5. Calcium 

Useful lines of Ca ii are found in the ESPaDOnS and HARPS 
spectra, in the near ultraviolet at 3933 A and in the ESPaDOnS 
infrared spectra at 8498, 8542, and 8662 A. The infrared lines 
are all blended with Paschen lines, and so the abundance of Ca 
was determined using only the /13933 line in the two ESPaDOnS 
spectra. The abundance of Ca appears to be approximately solar 
A weak line of Ca ii nearly coincides with a strong unidentified 
feature at 4799.9 A, but the computed strength of this line using 
the abundance derived from /13933 is negligible (see Figure[3]). 



5.6. Titanium 

Lines of Ti occur throughout the blue spectral region available 
in almost all the spectra. We have modelled lines at 4468, 4488, 
4518 A, and also Unes at 4798 and 4805 A. The two lines in the 
longer wavelength window are both well fit at all phases (see 
Figure O, but in the shorter wavelength window the computed 
line is somewhat too strong (and not quite the correct shape) 
for the strong line at 4468 A and somewhat too weak for the 
weak line at 45 1 8 A, which might be blended with another line 
not included in the line list. Nevertheless, synthesis of the two 
spectral regions yield very similar mean abundances, showing 
that Ti is overabundant with respect to the solar abundance by 
approximately 2 dex. 

5.7. Chromium 

The abundance of Cr has been derived from modelling lines of 
multiplets 39 and 44 at 4539, 4555, 4558 and 4565 A, and lines 
of multiplet 30 at 4812 and 4824 A. All lines are fit reasonably 
well with a uniform abundance. The various Zeeman split line 
profiles are approximately reproduced (see Figure [3]l. An over- 
abundance of almost 2 dex with respect to the solar abundance 
is found. This result is roughly consistent with the variations 
of abu ndance with T^ff found in Ap stars bv iRvabchikova et al.l 
(12004 . 

5.8. Manganese 

The strongest visible lines of Mn are in the region below the 
Balmer jump, between 3430 and 3500 A. Since these lines are 
not included in any of the available spectra, we modelled several 
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lines of Mn ii between 6122 and 6130 A. These high-excitation 
lines are not detected, so they only provide an upper limit of 
about 0.6 dex above solar abundance. It is interesting that al- 
though the even-Z iron peak elements Ti, Cr and Fe are over- 
abundant by 1-2 dex, this does not seem to be the case for Mn. 

5.9. Iron 

Iron has many usefully strong lines scattered through the spec- 
trum. We have fit lines in three windows, 4488 - 4542, 4820 - 
4830, and 5026 - 5037 A. All the numerous Fe ii lines in these 
regions are reasonably well fit with an abundance about 1.5 dex 
above the solar abundance. The weaker lines seen in Figure |4] 
suggest that the abundance of Fe might be even a little larger 
than this value; in any case, they show the tendency, also found 
for Ti, that computed weak lines to be weaker than the observed 
lines. This result is roughly with the variation of Fe ab undance 
with Teff found in Ap stars bv lRvabchikova et al.l (|2004|) . 



5.10. Praseodymium 

Four useful lines of Pr iii are found in the region 6098 - 6196 A. 
All are fit reasonably well with a mean abundance almost 5 dex 
larger than the solar value. 



5.11. Neodymium 

Nd m has lines scattered through the visible window. We have 
modelled Unes at 4796, 4822 (Figure|3]l, and 5050 A(Figure|4li, 
and find a mean abundance about 4.3 dex larger than the solar 
abundance. 



6. Discussion & conclusions 

This paper is the first in what is expected to be a series of reports 
describing modelling of stars from the survey of magnetic Ap 
stars in open clusters of known age. The aim of this project is to 
establish a preliminary model of the magnetic field structure, to 
estimate the chemical abundances of a number of common el- 
ements, and to extract whatever information is available in our 
data about the surface homogeneity or inhomogeneity of those 
elements by means of simple, rather schematic models whose 
parameters are derived by fitting synthesised spectra to observa- 
tions. 

HD 318107 is a relatively hot magnetic Ap star, with Teff - 
11 800 K. It has log L/Lq = 1.92, and mass M/Mq = 2.95. The 
rotation period of the star is 9.7088 d, and it has vsin/ = 7 km 
s"'. From its membership in NGC 6405 we know that its age 
is log t = 7.80 and that it has completed about 17% of its main 
sequence lifetime. The star has a strong magnetic field, of order 
15 kG globally. 

We have used pre viously pub lished (and modelled) magnetic 
field measurements of Mathvs et al. ( 1997), Manfroid & Mathvs 
(l2000l) . lLandstreet & Mathvsl (12000). and Bagn ulo et al.. (.2006) . 
together with measurements of our own, to improve the rotation 
period of the star to f = 9.7088 + 0.0007 d, which makes it 
possible to phase recent magnetic data and spectroscopy accu- 
rately with previous measurements, and thus to use the full set 
of magnetic and spectral data available to us to model the star. 

Our model of the magnetic field is a simple, low-order ax- 
isymmetric multipole expansion which retains the global topol- 
ogy of a dipolar field structure. This model is used in the frame- 
work of the oblique rigid rotator model to describe the observed 



variations of (Z?,) and (B) with rotational phase. Fitting the pre- 
dicted (B) and (Z? ) variations approximately to the observed 
variations allows us to determine the parameters of the model 
magnetic field as well as the inclination / of the rotation axis 
and the obliquity /? of the magnetic field to the rotation axis. 
Because the variations of the field moments are rather complex, 
and in fact poorly determined for {BJ, the adopted model is only 
a rather coarse first approximation to the actual magnetic field 
geometry of HD 318107. However, this model is appropriately 
simple considering the limited information available about the 
star, in particular the small number of (B,) measurements and 
high-dispersion polarimetric spectra available. 

The qualitative form of the field structure is clear from the 
nature of the field data shown in Figure[T] The fact that (B.) does 
not reverse sign indicates that the line of sight to the star is con- 
fined to one magnetic hemisphere, which in turn means that the 
sum i+/3 < 90°. The fact that the minimum value of is con- 
siderably smaller than the maximum value shows that the line of 
sight does approach the magnetic equator, while the large ratio 
of the maximum value of \{B,)\ relative to the value of (B) indi- 
cates that the line of sight does approach the magnetic pole. The 
actual parameter values adopted are discussed in Sec 4. 1 . The re- 
sulting magnetic model does not reproduce particularly well the 
detailed shape of the variations of (B ), but it works surprisingly 
well as a description of the Zeeman broadening (the dominant 
source of broadening) of the spectral lines we analysed. 

With only a dozen I spectra scattered through the rotation 
cycle (but including spectra obtained near times of maximum 
and minimum {B,}; i.e. when the line of sight is near the mag- 
netic pole and near the magnetic equator), it is not possible to do 
detailed mapping. However, a preliminary reconnaissance of the 
abundances of several elements, and a first estimation of their 
variation over the surface, is possible. We have considered a 
model that uses three co-axial rings, one around each magnetic 
pole and one around the equator, but there is no strong evidence 
for strong abundance variations on the global scale. The model 
finally adopted assumes uniform abundance over the stellar sur- 
face. This model provides the simplest possible quantitative de- 
scription of abundance variations on the star. 

Assuming uniform abundances over the stellar surface, we 
find, as expected, that most elements have distinctly non-solar 
abundances. The abundance of O is slightly below solar, while 
Mg appears to be significantly lower than the solar abundances. 
Ca has nearly solar abundance. Upper limits to the abundances 
of He and Mn are possible, the former unambiguously classify- 
ing HD 318107 has a He-weak star and the latter suggesting an 
abundance no more than 0.6 dex above the solar value. The other 
studied elements (Si, Ti, Cr, Fe, Pr, and Nd) are all present with 
higher abundance than in the Sun. For all of the elements stud- 
ied, the abundances measured over the part of the surface visi- 
ble in our spectra do not appear to vary strongly with magnetic 
co-latitude. However, for some elements a fit to the five spectra 
modelled is significantly better using the three-ring model than 
with a uniform abundance; we interpret this result as reflecting 
substantial small-scale abundance variations, as already hinted 
by the very different (B,) values measured in the ESPaDOnS 
spectra using lines of different elements. 

HD 318107 appears to be a star which may merit further 
study. With its extremely large and apparently rather complex 
field, and large (and apparently fairly non-uniform) abundance 
anomalies, this may be a good candidate for more detailed map- 
ping, which would require a major observing campaign to obtain 
phase coverage in all the Stokes components. The main diffi- 
culties of such a study are the lack of much Doppler broaden- 
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ing of the lines to provide spatial resolution in azimuth, and the 
smaU value of i which limits the line of sight to one magnetic 
hemisphere. But even before this is done, obtaining a few more 
spectropolarimetric V observations should help to clarify the na- 
ture of the variation of {B^} with phase, enable more detailed 
modelling of the field geometry, and help us to understand the 
discrepancies between the values of (B^) measured with Balmer 
Unes and with various metals. 
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